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Abstract 



In the Minimal Supersymmetric Standard Model with bilinear R-Parity 
violation, terms that violate R-Parity and lepton number are introduced in 
the superpotential, and sneutrino vacuum expectation values are induced. As 
a result, neutrino masses and mixing angles are generated via a low energy 
see-saw mechanism. We show that this model embedded into an anomaly 
mediated supersymmetry breaking scenario is testable at a linear collider 
using charged Higgs boson production in association with a stau. This is 
possible in regions of parameter space where the charged Higgs and stau 
have similar mass, producing an enhancement of the charged scalar mixing 
angles. We show that the bilinear parameter and the sneutrino vev can be 
determined from charged scalar observables, and estimate the precision of this 
determination. 
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I. INTRODUCTION 



Over the last three decades experimental evidence has confirmed the gauge structure 
of the Standard Model (SM) with very accurate measurements. Nevertheless, the picture 
is still incomplete since the Higgs mechanism has not been established yet experimentally. 
This mechanism lies in the center of the mass generation problem, giving mass to the gauge 
bosons as well as the quarks and leptons. Despite this success, it is clear that the SM should 
be extended. Theoretically, the SM does not have an answer to the gauge hierarchy problem, 
nor to the stability of the Higgs mass under quantum corrections. Supersymmetry is one of 
the most popular extensions of the SM that addresses satisfactorily these problems. 

Even if supersymmetry is not a symmetry chosen by nature, there is a generalized feeling 
in the community that important discoveries are going to be available soon after the com- 
pletion of the new generation of colliders, starting with the Large Hadron Collider (LHC) at 
CERN, or maybe even at the already running Tevatron at Fermilab. A Linear Collider (LC) 
is crucial in order to study the new phenomena and its relation to physics at even higher 
energy scales. This relation is more than a simple complementarity with the LHC, since the 
ideal is for them to run simultaneously in such a way that the discoveries in one machine 
may influence the running parameters of the other [1]. 

In parallel to the physics we can learn form colliders goes neutrino physics. Today, 
neutrino physics has become one of the most exiting areas in particle physics, with exper- 
imental results that indicate neutrinos have a mass and oscillate. Atmospheric neutrino 
data indicates a — v T mixing with an angle 0.3 < sin 2 # atm < 0.7 and a mass splitting of 

0.03 < \j Aml tm < 0.07 eV. Solar neutrino data favors a LMA solution with the best fit given 

by tan 2 9 sol = 0.44 and y / Am^ = 0.008 eV [2]. There is little doubt now that neutrinos 
have mass and this is the first experimental evidence that the SM must be modified. 

Supersymmetry with Bilinear R-Parity Violation (BRpV) [3-8] is a model in which neu- 
trino masses and mixing angles are generated by the presence of bilinear terms in the super- 
potential which violate R-Parity and lepton number. At the same time, sneutrino vacuum 
expectation values are induced. In this model, neutralinos mix with neutrinos, generating 
at tree level a mass for one of the neutrinos, while the other two acquire a mass at one loop. 

One of the nicest features of supersymmetric models which violates R-Parity and lepton 
numbers through BRpV is that neutrino physics is closely related to high energy physics 
in such a way that neutrino properties can be tested at future colliders, in particular at 
a LC where it is possible to make precision measurements on different observables. It is 
already understood that neutrino mixing angles are related to ratios of branching ratios 
of the neutralino, which being the Lightest Supersymmetric Particle (LSP), will have non- 
suppressed R-Parity violating branching ratios [4,9]. 

In BRpV there is also mixing in the scalar sector, in particular, charged Higgs fields mix 
with charged sleptons. As a result, in electron positron collisions it is possible to produce a 
charged Higgs boson in association with a stau, for example. This kind of R-Parity violating 
processes are not present in trilinear R-Parity violating models and, therefore, would be a 
signature of BRpV. In this paper we study charged Higgs production in association with a 
stau and its relation to neutrino physics in a model with Anomaly Mediated Supersymmetry 
Breaking (AMSB) [10-12], where an enhancement of charged scalar mixing angles can occur 
due to mass degeneracy between the scalars [13]. We study how parameters relevant to solar 
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neutrino physics, like the parameters in the superpotential, can be extracted from collider 
observables in the charged scalar sector. 



II. BILINEAR R-PARITY VIOLATION AND NEUTRINO PHYSICS 

In the BRpV-MSSM model, explicit lepton number and R-Parity violating terms are 
added to the MSSM superpotential 

W = Wmssm + eiliHu (2.1) 

where the three q parameters have unit of mass. In addition, corresponding soft terms add 
included in the lagrangian 

£soft = £^,ff M + BieiLiHu (2.2) 

with Bi being the bilinear soft terms associated to Cj. These terms induce sneutrino vacuum 
expectation values Vi which contribute to the gauge boson masses. Our notation for the 
fields that acquire a vev is 

Hd J»}\, H J»i\ z t J% (, 3) 



where 



H^^ + v. + i^], ^eI^ + ^I L° = -L[v« + v i + ii>!]. (2.4) 
The tree level scalar potential contains the following linear terms 

^linear = t d (y d + K°u + h V \ + h V 2 + h V Z , I 2 - 5 ) 

where the different t° are the tadpoles at tree level. They are given by 

t°d = ( m H d + V 2 ) v d + V d D - h(Bv u + VjSj) 

€ = -B(iv d + (m 2 Hu + /i 2 )w u - v u D + VjBjSj + v u e 2 (2.6) 
t ■ = ViD + ti ( - fiv d + v u Bi + Vj€j^ + ViM^ 

where we have defined D = |(g 2 + g' 2 ){vf + v\ + v\ + v d — v 2 .), e 2 = e\ + e 2 , + e\, and there 
is a sum over j but not over i. The five tree level tadpoles i° are equal to zero at the 
minimum of the tree level potential, and from there one can determine the five tree level 
vacuum expectation values. 

As a consequence of the e terms and the sneutrino vevs R-Parity and lepton number are 
violated, and R-Parity even fields mix with R-Parity odd fields. The most conspicuous of 
these is the mixing between neutralinos and neutrinos. This is important because due to a 
low energy see-saw mechanism, the neutrinos acquire mass. 



M l9 2 + M 2 g' 2 

m eff = 



4det(M x o) 



A 2 A e A M A e A r 
A e A M A 2 A M A T 
A e A T A M A r A 2 



(2.7) 



where M. x o is the 4x4 neutralino mass matrix, and we define 



Aj = fiVi + 6iV d . (2.8) 

These important parameters define the neutrino physics at tree level. The effective neutrino 
mass matrix in eq. (2.7) has only one eigenvalue different from zero, and it is given by 
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A — ; — Mio 2 + M 2 q- , r 

A '"-^'"-° Lam ") |A| ■ (2 ' 9) 

Since the lowest two eigenvalues are massless and degenerate, there is no meaningful solar 
angle at tree level. On the other hand, the near maximal atmospheric angle is 

tan #23 = A^/A T « 1 (2.10) 
and the Chooz constraint sin 2 9 13 < 0.045 [15] can be satisfied taking 

|tan0i 3 | = (Ael/^ + A 2 < 1 . (2.11) 

Information on the solar neutrino mass and angle can be obtained only when one-loop 
corrections are added to the neutrino mass matrix. As an example, the bottom/sbottom 
loops give a simple contribution which is approximately given by 

^Am sol ~ m U2 ~ w j^— In ^ . (2.12) 

In supergravity models, this is the most important loop, followed by charged Higgs loops. It 
is clear from this formula that the e parameters have an important effect in solar neutrino 
physics, as opposed to the A parameters which are important for atmospheric neutrino 
physics. 

When details of the neutrino physics are not the main issue, it has been proven very 
useful to work in the approximation where BRpV is introduced only in one generation, say 
the stau. In this case, the atmospheric scale still is given by eq. (2.9). In the rest of this 
article we will follow this approach. 



III. BENCHMARK FOR BRPV-AMSB 

In order to study the effects of BRpV in an AMSB scenario, we have chosen a case study, 
or benchmark, in which we find non negligible charged Higgs production in association with 
a stau. The parameters which define our model are 

M 3/2 = 30 TeV 

m = 200 GeV 
tan (3 = 15 

H<0 (3.1) 
e 3 = 1 GeV 
m v = 0.1 eV 
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The value of /i is fixed by imposing the correct electroweak symmetry breaking to be /x = 
—466 GeV. The neutrino mass fix the sneutrino vacuum expectation value to v 3 = 0.035 
GeV. The values of B and B 3 are also determined by the tadpole equations. 

Regarding the spectrum of this model, the LSP is the lightest neutralino, with a mass 
94.39 GeV, followed closely by the lightest chargino with a mass = 94.43 GeV. In 
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the neutral scalar sector we have a CP-odd Higgs with ttia = 171 GeV, and a tau-sneutrino 
with my r = 142 GeV. It should be stressed that in BRpV models the CP-odd Higgs boson 
mixes with the sneutrinos. We here call the CP-odd Higgs to the eigenstate with largest 
component to the original Higgs fields. 

In this article we are interested in the charged scalar sector. The charged Higgs fields 
and the left and right staus mix to form a set of four charged scalars, one of them being 
the charged Goldstone boson. Among the other three eigenstates we call charged Higgs to 
the one with largest component to the original charged Higgs fields. The charged scalar 
spectrum is 



m~+ 
m~+ 



188 GeV 
115 GeV 
190 GeV 



(3.2) 



The 4x4 charged scalar mass matrix is diagonalized by a 4 x 4 rotation matrix, which for 
our benchmark is given by 



R 



-0.067 
-0.80 
-0.014 
0.60 



0.998 
-0.05 
-0.001 
0.04 



-0.0001 
0.45 
0.65 
0.62 




0.40 
-0.76 
0.51 



(3.3) 



where the columns correspond to the fields H£, f£, f#) and the rows to 

(G + , H + , ?i, f 2 + ). From this rotation matrix we learn that the Goldstone boson has no 
right handed stau component and very little left handed stau component, as it should be. 
We also see that the light stau has almost no component to the Higgs fields, i.e., it is almost 
pure stau, and that the charged Higgs has an important component of stau. 



IV. CHARGED HIGGS/SLEPTON SECTOR 



As we mentioned, in BRpV charged Higgs bosons and charged sleptons mix forming, in 
the general three generations case, a 8 x 8 mass matrix, and in the simplified case of BRpV 
only in one generation, a 4 x 4 mass matrix. The relevant mass terms in the scalar potential 
are 



quadratic 



H u, H d 



1 L > ' R 



M 2 



(4.1) 



In the R- Parity conserving limit (MSSM), the mass matrix is diagonal in 2 x 2 blocks, and 
the charged Higgs sector is decoupled from the charged slepton sector. We write the mass 
matrix in the following form 
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M 2 



M| ( ± 0) + M 2 ^ 



(4.2) 



motivated by the fact that BRpV terms are small. The MSSM part is 

2(0) 9 



M 



2(0) 
5± 



2(0) 






2(0) 









" 










M\ R 


Ml R 


Mr 3 -1 



(4.3) 



where the diagonal slepton mass entries are given by 



Ml 3 
Ml 
M 2 LR 



Ml-l(g 2 -g' 2 M 
mf\A T - //tan/3) 



2(0) 



2(0) 

T 1 



(4.4) 



and is the tau lepton mass calculated in the MSSM approximation. Note that since in 
our model the tau mixes with the charginos, the physical tau mass has a more complicated 
dependence on the parameters of the model [14]. We write the BRpV contributions to the 
charged scalar mass matrix as 



M 2 « 



' Am% 

ti u 





X U L 


X U R 





Am 2 Hd 


XdL 


XdR 


X U L 


XdL 


Am 2 L3 







XdR 





Am| 3 



(4.5) 



with the diagonal contributions given by: 



Am 2 



l, fo S3 _ 1^2 2 , l/ ) 2 ? .2 



Am 2 H =4iml 

a 8 



v 3 c _|_ 12 2 



(4.6) 



Ami 



±fl 2 7 ,2 

g9z V 3 ) 



lh 2 v 2 4- V 2 ?; 2 



where m~ 



m 



(0)2 



I„2„,2 
8' 



g\v 3 and (/| = g 2 + g' 2 . The quantity m- - 1 is the sneutrino mass in 



the e 3 — > 0,f 3 — > limit [13]. The off-diagonal BRpV terms of in the charged scalar matrix 
in eq. (4.5) are 



X u l = 4# v d v 3 - fie 3 - ^h 2 v d v 3 , 
X uR = -^=h T (A T v 3 + e 3 v u ) , 

X d L = % 

Vd Sp 

X d R = —^h T (pv 3 + e 3 v d ) ■ 



C P , 1 2 

S/3 



(4.7) 



The complete charged scalar mass matrix in eq. (4.2) has a zero eigenvalue corresponding 
to the charged Goldstone boson G ± . This eigenvalue can be isolated with the rotation 
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Rg = 



where we have defined the factor 





S/3T 


Vd 


sp 


c (3 





-c 2 B r 

Vd P 


V d S P C P T 


r 












(4.8) 



r = 



(4.9) 



In what follows we present some approximated formulas for masses and mixing angles 
between charged Higgs bosons and charged scalars. Our aim is to find second order correc- 
tions to the masses and first order mixing angles. In order to do that, it is important to 
perform first the rotation defined in eq. (4.8), otherwise, approximating over the unrotated 
4x4 matrix in eq. (4.3) would introduce fictitious corrections to the zero Goldstone boson 
mass. In addition we rotate the stau sector by an angle Of : 



after which the zeroth order mass matrix becomes diagonal: 



"1 








" 





1 














Cf 


Sf 


.0 





-Sf 


Cf. 



(4.10) 





m 



2(0) 
H± 







Am 2 H± 





Xh?i 



2(0) 



m~y + Am? 



Tl 



Xhti 

X Hf2 

The corrections to the diagonal elements up to second order are 

™2 ,2 a ^,2 i J2 a ~~2 





XHf2 





m 



2(0) 
T2 



+ Am? 



T2 J 



(4.11) 



Am H± = SpAm Hu + c p Am Hi , 
Am 2 fi = 4Am| 3 + s?Am^ 3 + 4c|c f (c f M\ 3 + s f M 2 LR 



d 

Ami = 4*™ 2 L3 + c 2 f Am 2 Rs + $c 2 s f (sfM 2 3 - CfM 2 LR 

d 



(4.12) 



and the off diagonal corrections up to first order are 

Xhti = c f { s pX u L + CpX dL ) + Sf (spX uR + CfsX dR ) , 
Xht 2 = ~ s f [SpXuL + CpXdL) + Cf (spXuR + CpXdR) . 



(4.13) 



The last mixings X Hfl and X Hf2 are rotated away with the aid of two small rotations 



Ri 



"1 








0" 





Cl 


Si 








-si 


Cl 














1 



R 2 





c 2 

s 2 






1 







S2 



c 2 



(4.14) 



with rotation angles given by 
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Si = 



X 

2(0) 



Hfi 



2(0) ' 



S2 = 



X 

2(0) 



Hf 2 



m. 



2(0) 
T2 



(4.15) 



These are the mixing angles that are potentially large when the Higgs and stau masses are 
nearly degenerate. The final rotation that diagonalizes the charged Higgs/stau mass matrix 
is given by 



R — RiR 2 Rf Rc 



(4.16) 



It is clear that these approximations will fail in the case that two eigenvalues are nearly 
degenerate. We call 



sin 



2q+ 



R-23 ~l" -R- 2 



24 



(4.17) 



which indicates how much of stau fields has the charged Higgs. In Fig. 1 we compare the 
approximated formula with the exact numerical value of sin# + . The relative error is plotted 
against tan/3 for fixed values of M 3 / 2 = 30 TeV, m = 200 GeV, e 3 = 1 GeV, m v — 0.1 eV, 
and /i < 0. The error stays within 5% except for the points of near degeneracy. Close to 
tan/3 = 15 the charged Higgs and heavy stau masses are similar, and for tan/3 pa 17 the 
charged Higgs and light stau masses are similar. 

The second order corrections to the Higgs and stau masses are 



Am 2 



Am 2 fi 



-_: \ 1 1 1 

Tl 

2 _ \"2 



Arh 2 H ± + SiX Hfl + s 2 X Hf2 , 



Am 2 - S!X Hfl , 



(4.18) 



Am f2 = Am f2 



and since they are of second order in R-Parity violating parameters, they are small, although 
they may be enhanced if there is near degeneracy. 

V. CHARGED HIGGS AND STAU ASSOCIATED PRODUCTION 

A couple of charged scalars are produced in electron-positron annihilation via the in- 
terchange in the s-channel of a photon and a Z boson. The total production cross section 
is 



a(e + e- - S+ffr) 



A 3/2 

247TS 



9e*9v D x \ij , 9 2 (9v 



n 2„e 



-Sii - 



2c 



w 



-PzSijXzss + 



p2 \ij2 



(5.1) 



with gy — | — 2sfy, g e A = |, and = s/(s — m 2 z ). In the case of mixed production, only 
the Z boson contributes, with a strength determined by the ZS + S~ coupling, X l zss- m the 
unrotated basis (H u , H d , f L , f R ) these couplings are 



a: 



zss 





—C-2W 








" 


g 





-C2W 








2cw 








C-2W 
















9 s 2 



(5.2) 
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Since ~ 1/4 the couplings are approximately proportional to \' Z ss ~ diag(l, 1, 1 — 1). 
The rotated couplings are of course Xzss = RA' Z55 R T . 

For the chosen benchmark, the non negligible production cross sections are given by, 



o-(e + e~ 




= 22.9 fb 


a(e + e~ - 


-> if+if-) 


= 22.2 fb 


a(e + e~ 




= 21.0 fb 


<r(e + e~ 




= 3.44 fb 


a(e + e~ 


- 


= 1.81 fb 






= 0.79 fb 



(5.3) 



where the last two violate R-Parity. With a projected integrated luminosity of L = 500 fb -1 , 
we expect plenty of events where charged Higgs are produced in association with staus. 
These cross sections are governed by the couplings A^ss, which for our benchmark are 
approximately 



A 



zss 





"-0.5 
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-0.24 


-0.39 


0.2 


2c w 





-0.39 


0.08 


-0.3 







0.2 


-0.3 


-0.34 



(5.4) 



written in the base (G + , H + , ff, f£)- The first thing to notice is that the Z couples to a pair 
of light staus with a strength 3 to 4 times smaller than to the other two scalars. The reason 
is that the light stau is almost pure stau and has nearly equal left and right component, and 
the minus sign in \' zss induces a cancellation. For this reason its cross section is comparable 
to the other two pair productions even though there is more phase space for light staus. The 
same type of cancellation does not happen to the heavy stau because this eigenstate has 
also a large component of Higgs fields. 

The second thing worth noticing is that despite the fact the light stau is almost pure 
stau, as indicated by eq. (3.3), its mixed production together with a charged Higgs is larger 
than the mixed production of a heavy stau together with a charged Higgs, as shown in 
eq. (5.3). The explanation is simple: the charged Higgs has a large component of stau. 

In Fig. 2 we plot as a function of tan f3 the expected number of events in which a charged 
Higgs is produced in association with a stau, where we have summed over both possible 
signs of the scalar electric charge. The cross section has a maximum near tan/3 = 15 where 
the charged Higgs mass is very similar to the mass of the heavy stau. The charged Higgs 
associated production with a light stau also has a maximum at tan (5 ~ 15 because at this 
point the charged Higgs has a large stau component, as explained before. As we increase 
tan/5 the charged Higgs mass decreases and eventually, near tan/9 16.7, becomes similar 
to the light stau mass, increasing again the cross section. 

In order to have an idea of the degree of degeneracy needed for large cross sections, we 
plot in Fig. 3 a similar graph but this time as a function of the mass splitting — mf 2 . 
A mass difference of 10 GeV gives of the order of 50 events, and increasing exponentially 
with decreasing mass difference. 

The dependence of these two cross sections on the scalar mass m can be seen in Fig. 4, 
where a maximum value is achieved at m = 190 GeV. Nevertheless, the maximum value 
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for the cross section in much smaller than in the previous two figures, where the control 
parameter is tan /3. 

To finish this section, in Fig. 5 we present the production cross section times luminosity 
of a charged Higgs in association with a stau as a function of the gravitino mass M 3 / 2 . A 
sharp maximum is observed for both cases at M 3 / 2 = 30 TeV and decreasing rapidly for 
smaller values. At higher values the cross sections start to rise again, but the curves end 
when too low values of scalar masses are reached. 



VI. EXTRACTING SUPERSYMMETRIC PARAMETERS FROM OBSERVABLES 

As it was mentioned before, ratios of branching ratios of neutralino decays are related to 
the parameters Aj, which are in turn related to the atmospheric mass and angle. In this way, 
measurements on neutralino decays may give information on these parameters. Nevertheless, 
it is not easy to extract from neutralino physics information on the 6j parameters or on the 
sneutrino vevs V{. 

The situation is different in the scalar sector, and in particular, in the charged scalar 
sector. As we have seen in chapter IV, mixing angles in the Higgs stau sector depend directly 
upon the BRpV parameter e 3 and the sneutrino vacuum expectation value v 3 . In addition, 
near degeneracy between the charged Higgs and the stau found in AMSB enhances the 
associated production cross section, making this R-Parity violating process observable. 

Based on this idea, we study the possibility of extracting the fundamental parameters 
of the model, and especially e 3 and t> 3 , from hypothetical measurements of production cross 
sections and decay rates of charged scalars. We use a x 2 method with an input given by the 
AMSB model defined in eq. (3.1). We randomly generate models varying all the relevant 
parameters in the MSSM: 

• Higgs sector parameters /i, tan/3, m^. 

• Slepton soft mass parameters M|, M|. 

• Trilinear soft mass parameters A T . 

• Gaugino masses M 1; M 2 . 

• Neutrino sector parameters e 3 , m U3 . 

These parameters are free, in the sense that are not calculated with the AMSB boundary 
conditions. For the ith model we calculate a xf (normalized by the number of observables) 
given by 

where <7j is the observable calculated with the random parameters of the ith model, a is the 
observable calculated with parameters given by the input AMSB benchmark, and bo is a 
projected error in the measurement of the observable, estimated using only the statistical 
error N±\/~N. As observables we use the charged scalar masses in eq. (3.2) (with a projected 
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error of 1%), the production cross sections in eq. (5.3), and branching ratios of charged scalars 
decaying into charginos, neutralinos, and leptons. 

From a large sample of models we select 2000 of them with normalized x 2 < 2 and 
plot the ^-distributions in Figs. 6 and 7. The most interesting distributions for us are the 
corresponding to e 3 and V3 in Fig. 6. In this figure we see that a clear upper bound on 
these parameters can be set, in addition to a less clear lower bound. This is an important 
achievement considering the difficulties in extracting values for these parameters from neu- 
tralino decay measurements. From the rest of the ^-distributions we learn that a reasonable 
determination of tan/5, m^, M^, M^ 3 , /i, and M 2 can be made. No useful information can 
be obtained for the parameter A T , and only limited information for Mi (which we do not 
show) . 

In Fig. 8 we have regions in different planes of parameter space where normalized x 2 < 1 
(light grey, or green), which tells us about the error in the determination of the corresponding 
parameter. For comparison we also show the regions with x 2 < 2 (dark grey, or blue). From 
this figure we extract the output for each parameter from the x 2 analysis, the error in 
its determination, and compare them with the input values from our benchmark. This 
comparison is shown in Table I. The determination of all the parameters of the R-Parity 
conserving MSSM directly involved in the Higgs/slepton sector is very good, with a few 
percent of error. The fact that we can set an upper and lower bound on the R-Parity 
violating parameters e 3 and v 3 is good in itself, although with large errors. 

The enhancement of the charged scalar mixing angles due to near degeneracy will modify 
the prediction of the solar mixing angle and mass scale. This possibility will be studied in 
a further work where BRpV in three generations is included. 
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VII. CONCLUSIONS 



Supersymmetric models with bilinear R-Parity violation predict neutrino masses and 
mixing angles which agrees with experimental data on solar and atmospheric neutrinos. 
It has been shown that the supergravity version of this model is testable at colliders via 
neutralino decays, when this particle is the LSP, and that information on the parameters A 
can be obtained. In this article we show that BRpV embedded into an Anomaly Mediated 
Supersymmetry Breaking model can also be tested at colliders with processes not necessarily 
related to LSP decay. In fact, we show that in cases of near degeneracy between charged 
Higgs boson and staus, the charged scalar mixing is enhanced, as well as the charged Higgs 
production in association with staus. The end result is that it is possible to determine the 
values of the parameters €3 and v 3 from measurements of charged scalar masses, production 
cross sections, and decay rates. 

Note Added: When this article was being written we received a related work where neu- 
trino physics is probed at colliders with charged slepton decays in Supergravity when the 
slepton is the LSP [16]. The ideas presented in their article are in agreement with ours. 
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TABLES 

TABLE I. Parameter determination from charged Higgs and stau associated production and 
decay for a given case study (benchmark). All parameters are expressed in GeV, except for tan/3. 



parameter 


input 


output 


error 


percent 




1.0 


1.24 


0.96 


77 




0.035 


0.045 


0.035 


78 


tan/3 


15.0 


15.0 


0.3 


2 


A* 


-466 


-469 


22 


5 


Ml 3 


155.8 


155 


3 


2 


Mr 3 


144.5 


144 


3 


2 


m A 


171.5 


172 


3 


2 


M 2 


95.4 


95 


3 


3 
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FIG. 1. The relative error of the approximated formula for the angle sin0 + , which indicates 
how much component of stau fields has the charged Higgs, is plotted against tan/3. 
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FIG. 2. Charged Higgs production in association with a stau multiplied by an integrated lumi- 
nosity of 500 ftT 1 , plotted as a function of tan/3, for M 3/2 = 30 TeV, m = 200 GeV, e 3 = 1 GeV, 
m u = 0.1 eV, and (jl < 0. 
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FIG. 3. Charged Higgs production in association with a stau multiplied by an integrated lu- 
minosity of 500 fb~ , plotted as a function of the mass splitting — mf 2 , for M 3 / 2 = 30 TeV, 
m = 200 GeV, e 3 = 1 GeV, m v = 0.1 eV, and p<Q. 
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FIG. 4. Charged Higgs production in association with a stau multiplied by an integrated lu- 
minosity of 500 fb , plotted as a function of the scalar mass mo, for M 3 / 2 = 30 TeV, tan/3 = 15 
GeV, e 3 = 1 GeV, m u = 0.1 eV, and fj, < 0. 
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FIG. 5. Charged Higgs production in association with a stau multiplied by an integrated lumi- 
nosity of 500 fb , plotted as a function of M 3 / 2 , for tuq = 200 GeV, tan/3 = 15 GeV, €3 = 1 GeV, 
m v = 0.1 eV, and \i < 0. 
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FIG. 6. Normalized x 2 lower bound as a function of the parameters (a) e^, (b) v$, (c) tan/3, 
and (d) itla- The boundaries indicated by x 2 = 1 defines the expected error in the determination 
of these parameters. 
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FIG. 7. Normalized x 2 lower bound as a function of the parameters (a) Ml 3 , (b) Mr 3 , (c) /j,, 
and (d) The boundaries indicated by x 2 = 1 defines the expected error in the determination 
of these parameters. 
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